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The electronic properties of a graphene sheet with attached hydrogen atoms is studied using a 
modified Fahcov-Kimball model on the honeycomb lattice. It is shown that in the ground state 
this system separates into two phases: fully hydrogenated graphene (graphane) and hydrogen-free 
graphene. The graphene-graphane boundary acquires a positive interface tension. Therefore, the 
graphene-graphane interface becomes a straight line, slightly rippled by thermal fluctuations. A 
smooth interface may be useful for the fabrication of mesoscopic graphene-based devices. 
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I. INTRODUCTION 

Creating a sample with flat edges is a significant chal- 
lenge for producing graphene mesoscopic devicesi. One 
possibility is to break a graphene sheet into fragments 
with sharp edgesi^ Another alternative involves the use of 
graphane. Graphane^ is fully hydrogenated graphene; it 
is an insulator with a gap of several cV. With graphane, 
instead of physically cutting graphene, one can create 
graphene patches of required shapes inside a sheet of 
graphane by local de-hydrogenation. In such systems, 
low-energy electrons from graphene cannot penetrate 
the insulating graphane host. Therefore, the graphene- 
graphane interface serves as the effective edge of the 
graphene structure. Different arrangements of this type 
have been discussedi^ Thus, the issue of the graphene- 
graphane interface stability is important both for funda- 
mental and applied research. There are indications from 
numerical studies that such interface is stable^>i, and that 
the adsorbed hydrogens tend to cluster together^ This 
tendency may be explained in terms of phase separation 
into hydrogen-rich and hydrogen-free regions, which was 
established on the basis of the semi-phenomenological 
analysis of the electron-mediated interaction between hy- 
drogen adatoms in graphene 

The purpose of the present paper is twofold. First, 
we put forward a microscopic approach to the problem 
of phase separation in graphene-graphane systems. To 
demonstrate the phase separation Ref. assumed a spe- 
cific type of interaction between the graphene electrons 
and the adatoms. For the hydrogen adatoms this as- 
sumption is supported experimentally and numerically. 
However, it remains unclear if the phase separation is 
a unique feature of the hydrogen on graphene, or other 
adsorbents will show the same feature. Avoiding phe- 
nomenological arguments we discuss the phase separa- 
tion within the framework of a modified Falicov-Kimball 
model with an infinite interaction constant between 'a 
hydrogen hole' and an s-electron on the hydrogen atom. 
The advantage of such an approach is its generality: the 



phase separation is a known property of a ground state 
of Falicov-Kimball-like modelsiSrii robust against vari- 
ation of microscopic details. To estimate the character- 
istic energies of the phase-separated state we apply the 
Hubbard-I approximation^^ To check the validity of this 
approximation, we also perform exact diagonalization of 
the model Hamiltonian in a finite cluster. 

The phase separation implies that the homogeneous 
state is either unstable or metastable. However, it is 
possible to imagine that, under suitable conditions, such 
phase may be stabilized for substantial amount of time. 
If the stabilization is indeed possible, the properties of 
the homogeneous phase can be investigated. Our calcu- 
lations show that the homogeneous phase is metallic, in 
agreement with the numerical results of Ref. 

Our second goal is to explore the connection be- 
tween phase separation and the stability of the graphene- 
graphane interface. We show that the graphene-graphane 
interface has a positive boundary tension. To stretch the 
interface with a positive interface tension by a unit length 
requires a finite amount of work. This amount is high for 
the system considered. Thus, the interface remains flat 
over substantial distances, which is a highly desired prop- 
erty, necessary for the creation of ballistic mesoscopic sys- 
tems. In other words, the interface is stable not only with 
respect to vacancy defects in small samples, as found in 
Refs.i and0, but also with respect to any conceivable 
defect. Our approach allows to obtain a qualitative es- 
timate of the interface tension and to assess the flatness 
of the interface at a given temperature. We estimate 
that at room temperature the graphene-graphane inter- 
face remains atomically smooth over distances of about 
10^ lattice constants. 

The paper is organized as follows. In Scc.|ll]we formu- 
late the model of the adatoms adsorbed on a graphene 
sample. This model is solved in Sec. IIIII within the 
Hubbard-I approximation. To check the accuracy of 
the Hubbard-I approximation the flnite-cluster numerical 
study is presented in Sec. |lVl In Sec. |V] we investigate 
the stability of the graphene-graphane interface, evalu- 
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ate the interface tension, and investigate its smoothness. 
The conclusions are given in Sec. IVII 



II. MODEL 



We use the model Hamiltonian for graphane: 



Ha = He-J2 h {PlS^<r + h.c.) + eaSlS,^ 



,(1) 

(2) 



where = (pf;,pfj), Sl = (sf^.^fj)' ^ the spin 
projection. The Hamiltonian ffs (^^a) corresponds to 
graphene (graphane). Below, label 'E' ('A') is used to 
denote quantities associated with graphEne (graph Ane). 
The Hamiltonian i/s is the usual graphene Hamiltonian 
corresponding to p^-electrons of carbon hopping between 
nearest carbon atoms arranged into the honeycomb lat- 
tice. For such lattice, the electron creation operators 
are arranged into a spinor Pj^, where i denotes the bi- 
atomic unit cell of the lattice. The spinor component 
labeled ('S') corresponds to a site on the A sublat- 
tice [B sublattice) . The hopping matrix in the spinor 
representation in momentum space is 







1 + 2 exp 



cos 



(3) 



The Hamiltonian Ha is a simplified model of graphane. 
It describes the pj-electrons of graphene hybridized with 
the s-electrons of hydrogen, attached to each carbon 
atom. Other bands are disregarded. The carbon- 
hydrogen hybridization constant = 5.8 eV exceeds the 
carbon-carbon hopping amplitude tp = 2.7 eV and the 
relative energy of the hydrogen s-orbital eh = 0.4 eVi^ 
In k-space the Hamiltonian Ha can be expressed as 



Ha 



f 7k ^o^o 



(4) 



where ctq is the 2x2 unity matrix. Here the upper left 
2x2 corner corresponds to the carbons atoms, lower right 
2x2 block corresponds to the hydrogens, the remaining 
blocks describe the C-H hopping. 

The matrix for Ha is easy to diagonalize. As a result 
we obtain four graphane bands: 



^ (±|tk|± 



l^kl 



1,2,3,4. (5) 



In this formula Ejj is neglected for it is small. 

Although only four bands in graphane are considered 
in our model Hamiltonian, Eq. ([l}, it captures the main 



features of graphane: at half filling, Eq. ([T]) describes 
an insulator with a gap located at the V point. The 
value of the gap Eg (for the parameters written above 



Eg = 6.0 eV) is found to be consistent with Ref. Il6lll 
Note, however, that there is no consensus about the ex- 
act values of the graphane model parameters. But high 
precision is not important for the qualitative results ob- 
tained below. 

The Hamiltonian Ha is valid when the numbers of hy- 
drogen and carbon atoms are equal. If at some site hy- 
drogen is absent, then the hydrogen s-orbital is not avail- 
able for the electrons. This constraint may be enforced 
by introducing an infinitely-strong repulsion between the 
"hydrogen hole" and the electron on the s-orbital: 



Hea^Ha+uJ2sIn^''S,^, 



N.. 



hh 



(6) 
(7) 



where U — > +oo, and n^jg^ are the numbers of hydrogen 
holes at site i. These numbers can randomly take the 
values or 1 with mean value {n^^igj} = n^^, where n^^^ 
is the concentration of hydrogen holes per carbon atom. 
The Hamiltonian Hea is a version of the Falicov-Kimball 
model in which mobile p- and s-electrons interact with 
immobile "holes" whose concentration n^^ is fixed exter- 
nally. Thus, n^^ = 1 refers to graphene, n^^ = refers 
to graphane. Below, we will study partial hydrogenation: 
< n'^'^ < 1. 



III. CALCULATIONS 

An analogy between H^a, Eq. ([5]), and the Falicov- 
Kimball model is very useful for our purposes since the 
latter model experiences phase separation in a broad 
range of parameters j^Sril The reasons for the existence 
of phase separation here can be understood with the help 
of simple arguments. Since exceeds tp, let us study the 
limit 



to ^ tp. 



(8) 



We now introduce the electron operators a, b diagonaliz- 
ing those terms of H^a which do not involve the carbon- 
carbon hopping: 



Hea -He = J2 " + 

a 



(9) 



where 



V2 



(l-n'^)+n^Hg, 

(i-n^') + <V. 



(10) 
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We omit the sublattice and spin labels since the expres- 
sions are the same for any A, B, and a. The index a labels 
individual carbon atoms (i, j label unit cells). In Eq. ^ 
we neglect the term proportional to £h since eh <SC tp 
<o)- It follows from Eq. ([9|) that the on-site energy of 
the fermions a is much higher than the on-site energy of 
h for any n'^^, since Iq^U ^ tp. Thus, to lowest order in 
tp/toy these states are empty, and can be neglected. In 
this approximation 

-| E ^"^^[1 + ^("'' + ^f) + (11) 

(12) 



{af3) 

7 V2 - 1 w 0.41 



where (...) denotes summation over the nearest neigh- 
bors. From this equation we see that to separate two hy- 
drogen holes sitting on neighboring sites one must spend 
an energy of the order oitp^"^ {h]jjfj) . This corresponds to 
the attraction between the hydrogen holes (and between 
the hydrogen atoms) as in the model used in Ref. This 
attraction induces the phase separation. The additional 
correlations between the adsorbed adatoms (e.g., due to 
bond reorganization in graphane), which our model ne- 
glects, may be incorporated as an effective short-range 
attraction between the hydrogens. The effect of this at- 
traction is obvious: it favors phase separation. 

Using 77ea , Eq- © , we can derive the equation of mo- 
tion for the single-electron Green's function in the (w,k) 
representation: 



{lo + n)Gpp + Ti^Gpp + toGsp — 1, 



{uj + ^1 + e}i)Gsp + toGp 



UF,, 



0. 



(13) 



Here /i is the chemical potential, Gpp_sp and F^p are the 
Fourier transforms of the time-ordered propagators 

Gpp{i-j,t) = -i(TP,„{t)Pl{Q)), 
Gsp{i-i.t) = -i(TSUt)P}M). (14) 
Fsp{i~j,t) = -i(TNl^'\t)S,,{t)PliO) 

The propagator F^p requires an additional equation of 
motion, which relates F^p with the propagator 



Fpp = -t{TN^\t)P,,{t)Pl{0)). 



(15) 



To truncate the infinite set of equations for the Green's 
functions, we apply the Hubbard-I approximation. It 
is a simple mean-field scheme suggested in the semi- 
nal papersi^. The applicability of Hubbard-I and re- 
lated approaches has been tested in many cases (see, e.g., 
Refs. imili). In the Hubbard-I approach, Fpp is ap- 
proximated by the product 



> 



Hubbard-I 

y 
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FIG. 1: (Color online) Electron energy E vs the concentration 
71^ of hydrogen adatoms calculated in the Hubbard-I approx- 
imation (green solid curve) and by exact diagonalization of 
10 X 10 unit cells cluster (red dashed curve). The negative 
curvature of E{n^) over the whole range of is an indica- 
tion of the instability of the system toward the macroscopic 
separation into phases with = and = 1. The Maxwell 
construction is shown by blue dot-dashed line. The inset 
shows the energy difference between homogeneous and phase 
separated states calculated in the Hubbard-I approximation 
(green solid curve) and by exact diagonalization of 10 x 10 
and 6x6 clusters (red dashed and blue dotted curves, respec- 
tively). The model parameters are: to = 5.8eV, tp = 2.7eV, 
Eh = 0.4 eV. For exact diagonalization, U = 400 eV. 



This closes the system of equations (fT3|) . whose solution 
may be written explicitly as 



Grj 



Gs 



to + ^i + en 



{uj + I^L + £h)(w + fi + Tk)- n^tl 



, (17) 



where = 1 — n''^ is the hydrogen concentration per 
carbon atom. These equations are obtained in the limit 
U — !■ oo. Similarly, the Green's function 



GU^-J,t)^-i{TS,,it)Sl{0)) 
is calculated 



[ 1 ~ toGsp 
uj + fi + en 



(18) 



(19) 



{N^'')Gp 



(16) 



In the limiting case = (n^ = 1), the Green's 
functions in Eqs. (fT7|) and (fT9|) coincide with the exact 
Green's functions corresponding to the Hamiltonian 
of graphene {Ha of graphane). 

When the Green's functions are known, the density of 
states, the electron concentration, and the energy can 
be calculated as a function of fi. Fixing the electron 
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FIG. 2: (Color online) Electronic dispersion e(k) of the un- 
stable homogeneous metallic phase for = 0.2. Four single- 
electron bands, found with the help of Hubbard- I^^'^° are plot- 
ted for different points of the Brillouin zone. The chemical 
potential fj, is marked by the horizontal green dashed line; 
fi was calculated self-consistently to ensure that the electron 
concentration is {1 + n^) per carbon atom. The gap between 
the conducting and the valence bands at = 0.2 is smaller 
than the graphane gap (n^ = 1). 



concentration (1 + n^) per carbon atom, we find /i = 
/i(n") and the energy E = E[fi{n^)] at T = 0. 

The Hubbard-I results are presented in Fig. [TJ The 
energy-versus-density curve has negative curvature for 
any n^. This indicates the instability of the homoge- 
neous phase toward the phase separation. The energy 
of the phase separated state can be found with the help 
of the Maxwell construction. In our case, it is simply a 
straight line connecting the energy of the pure graphcne 
at = and the energy of the fully hydrogenated 
graphane at = 1. This means that the separated 
phases are pure graphene and pure graphane. 

The single-electron band structure of the unstable 
mixed graphene-graphane phase is shown in Fig. [2] It 
is interesting to note that the homogeneous phase is a 
metal. This is consistent with numerical results for small 
clusters 1^ 



IV. NUMERICAL CALCULATIONS 

To check our analytic approach, we also perform ex- 
act diagonalization of the Hamiltonian ^ on a finite 
honeycomb cluster containing 10 x 10 unit cells (200 car- 
bon atoms). Periodic boundary conditions are used. For 
each n^, hydrogen atoms are randomly distributed on 
the cluster, and we calculate the system energy by aver- 
aging it over 1500 configurations. 

To check the reliability of the numerical results we in- 
vestigate their dependence on the number of sites in the 
cluster (A^sitcs) and the number of the configurations used 
for the averaging (iVconfig)- In Fig. |3] the averaged energy 
(normalized per site) 
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FIG. 3: (Color online) The lower (brown) curve represents 
the averaged energy (E), Eq. (|20[l. The upper (blue) curve 
represents the dispersion, Eq. (|21l) . The curves are plotted as 
functions of the number of the disorder realizations A'^conflg- 
Both curves demonstrate saturation for A'conflg ^ 750. 



and the normalized energy dispersion 



(21) 



config-^ sit 



(20) 



are plotted as functions of A'config- In Eq- (|20p index Q 
labels different disorder realizations, E\Q] is the energy 
for a given disorder realization <d. Both {E) and D-^ 
demonstrate saturation for TVconfig ^ 750. This suggest 
that iVconfig = 1500 we used in our numerical calculations 
is sufficient to obtain reliable results. 

In Fig. d] the same quantities are plotted versus iVgitcs- 

— 1/2 

The dispersion decays as A^gjtcs • This means that the rel- 
ative strength of the energy fluctuations decreases when 
the cluster size grows, and the energy experiences the self- 
averaging. The energy itself saturates for A'sitcs ^ 125. 
Therefore, our choice of A^sitos = 200 is adequate. In ad- 
dition, the ratio De/{E) sets the relative error for (E). 
For A'sites = 200 this error is a fraction of a percent. We 
conclude that our numerical calculations are reliable. 

The most important results are shown in Fig. [TJ ob- 
serve the negative curvature of the function E{n^). It 
implies that the system is unstable and phase separates 
in two phases: with = (graphcne) and with ~ 1 
(graphane) . Unless is close to or 1 , the energy gain 
due to the phase separation is of the order of lO'^ K, see 
the inset of Fig. [TJ Thus, even at room temperature we 
can safely use the results obtained at zero temperature. 

Further, the numerically evaluated energy is of the 
same order as the Hubbard-I energy: the magnitude of 
the Hubbard-I energy is approximately two times higher 
than the numerical estimate (see inset in Fig. [TJ). Thus, 
the qualitative consistency between the numerical calcu- 
lations and the Hubbard-I results provides firm support 
to the ffirdings of Sec. IIIII 
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f Clusters from 3x3 to 10x10 unit cells. 

Averaging over 1500 configurations at n'=hflN =0.5 
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FIG. 4: (Color online) Finite size effects. The solid (green) 
curve connecting skew crosses represents the averaged energy 
{£), Eq. It saturates for iVsitcs > 125. The solid (red) 

curve connecting filled circles represents the energy disper- 
sion, Eq. (|2H) . It decays as N~^^^^ . The decay of the dis- 
persion implies that for large samples the energy is a self- 
averaging quantity. The value of the dispersion may be used 
to evaluate the accuracy of the estimated value of the energy. 
For Nsitcs ~ 200 the error for (E) is a fraction of a percent. 



INTERFACE TENSION AND INTERFACE 
STABILITY 



In the phase-separated state there is a boundary be- 
tween graphene and graphane. The geometry of the sta- 
ble inhomogeneous state depends on the sign and the 
value of the interface tension ctq- If co < 0, then the inho- 
mogeneous phase breaks into small clusters to maximize 
the boundary length. In the case of ctq > 0, the inter- 
face tension acts to minimize the length of the graphene- 
graphane border. In the case of a long strip this bor- 
der is a straight line (if the concentration of the hydro- 
gen adatoms is not small). However, at finite temper- 
atures, even for positive ctoi small thermal fluctuations 
destroy the perfect smoothness of the boundary between 
the two phases. The difference in the lattice symmetry 
between graphene and graphane, at the level of the elec- 
tron model, manifests itself through the values of the or- 
bital overlaps. In the model considered here we make an 
approximation regarding the orbital overlaps: we assume 
that several of them are equal to zero. 

Further, we neglect the difference between lattice con- 
stants in graphene and graphane. The contribution of 
the electron-electron interaction to the interface tension 
is also disregarded (we briefly discuss the effect of the in- 
teraction below). In other words, the interface tension in 
our model arises only due to the electron motion through 
the graphene-graphanc boundary. 

These assumptions can be justified post factum: (i) 
from our model it follows that the binding energy be- 
tween a hydrogen atom and the graphene-graphanc in- 



terface is of the order of tp, which is consistent with the 
results presented in Ref. 0; (ii) we pointed out above that 
the value of the graphane gap in our simplified model 
Eq. (dJ turns out to be consistent with other studies; 
(iii) we found that for intermediate hydrogenation the 
stable homogeneous phase is a metal, in agreement with 
Ref. B 

We will now evaluate (Jq in the limit Eq. ([8]). In this 
approximation, electrons in graphane are localized on the 
C-H valence bonds [see Eq. ([9])] and their contribution to 
(To is small (this contribution is proportional to tp/to and 
en). In graphene, electrons are moving from one carbon 
atom to its nearest neighbors. However, the electrons 
from graphene cannot penetrate into graphane since they 
have to overcome the graphane gap, which, according 
to Eq. dll), is of the order of to when Eq. (|8]) holds. 
Thus, each carbon-carbon bond connecting an atom in 
graphene with an atom in graphane does not contribute 
to the graphene electron kinetic energy. This, in effect, is 
equivalent to an increase in the kinetic energy of the elec- 
trons in graphene. The longer the interface, the larger 
the number of "broken" bonds. Thus, 



ao 



(22) 



where £b is the kinetic energy for each carbon-carbon 
bond. The numerical coefficient 



_ / l/Va « 0.6 for zigzag 
2/3 « 0.7 for armchair 



(23) 



characterizes the linear concentration of the carbons on 
the interface. The kinetic energy per bond is equal to 



£b 



So\tk\ 
(27r)2 



(24) 



where the integration is performed over the first Brillouin 
zone. 5*0 — 3^/30.^/2 is the area of the graphene unit cell, 
the factor 2 corresponds to two spin projections, (1/3) 
enters since there are three bonds in a graphene unit 
cell, and ik is defined by Eq. After integration, we 
have Eb ~ 1.05 tp, and in our approximation, 

(To ^ 0.6 tp/ao- (25) 
A more accurate calculation (following Ref. [2l| ) provides 



(To ~ 0.2tp/ao ~ 0.6eV/ao. 



(26) 



In the calculations presented above the contribution of 
the electron-electron interaction to ao is disregarded. 
The detailed account of the interaction goes beyond the 
scope of the present study. Yet, we would like to of- 
fer two observations. First, the contribution due to the 
interaction is of the same order as ao. Indeed, the lat- 
ter originates mostly from the energy of C-C bond. The 
chemical energy of C-H bond is of the same order (few 
eV). Thus, there is no energy scale in the system which 
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would be able to generate an overwhelmingly large con- 
tribution to the interface tension. Second, the contri- 
bution due to the interaction increases the tension. To 
prove this, let us neglect the interaction in the bulk, as 
it is usually done for graphene, but retain the interaction 
term for the electrons near the graphene-graphane edge. 
This assumption mimics relative importance of the in- 
teractions for electrons in lower dimensions. It is known 
that the repulsive interaction gives positive contribution 
to the electron energy (see, e.g.. Sec. I, § 6 of Ref . [2^ and, 
consequently, to the interface tension. 

We neglect the effect of the temperature T on the phase 
separation since the characteristic energies of the problem 
are much higher than ksT for any realistic T. However, 
the temperature fluctuations could affect the smoothness 
of the graphene-graphane interface even under such con- 
ditions. Following Ref. [2^ we can express the average 
square fluctuation of the deviation u of the interface hav- 
ing a length L as 



ksTL 
2'Kao 



Thus, we obtain 



knT 



(27) 



(28) 



to the interface tension due to the difference between 
lattice constants in graphene and graphane. 



VI. CONCLUSION 



We mapped the model of hydrogen atoms adsorbed on 
graphene on a Falicov-Kimball-like model. We demon- 
strated that this system has a strong tendency to phase 
separate. The thermodynamically stable state is inho- 
mogeneous: all adatoms cluster together, forming two 
phases: hydrogen-saturated graphane and hydrogen-free 
graphene. The interface between these phases has fi- 
nite and positive interface tension, which means that the 
boundary is stable and flat (if the number of hydrogen 
adatoms is not small). The estimated value of the in- 
terface tension is high and, at room temperature, the 
interface remains atomically flat over distances of about 
10^ lattice constants. This result may be of interest for 
fabricating graphene mesoscopic devices with weak edge 
scattering. 
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